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We consider a general unitary operator acting on two 
qubits in a product state. We find the conditions such that 
the state of the qubits after the action is as entangled as pos- 
sible. We also consider the possibility of using ancilla qubits 
to increase the entanglement. 
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I. INTRODUCTION 

Entanglement is a quantum mechanical feature 
that can be employed for computational and communi- 
cation purposes. During the last few years a big effort 
has been done in order to create entanglement in several 
laboratories ||] . This entanglement can then be used for 
many fascinating things such as teleportation Q], quan- 
tum cryptography and quantum computation In 
some of these experiments the entanglement is produced 
by starting out from a product state of two systems (typ- 
ically qubits) and using some physical process that gives 
rise to an interaction between them. Thus, one of the 
relevant problems in this context is to find ways of gen- 
erating "as much entanglement as possible" for a given 
experimental set-up, i.e. a non-local interaction. 

The first steps to answer this problem have been given 
in Rcf. (7|-|^. In particular, given a non-local Hamilto- 
nian Diir et al have found the optimal way of generating 
entanglement. It consists of applying some fast local op- 
erations during the interaction processes in such a way 
that the rate at which entanglement increases is always 
maximal. In some situations, however, one cannot ap- 
ply fast local operations during the process, but rather 
a fixed quantum gate is given. In this work we find 
the states \4))a and \4')b for which the entanglement 
of Uab \ 4')a \ iP)b is maximal, where Uab is an arbitrary 
unitary operator. Thus, our results give a characteriza- 
tion of two-qubit gates in terms of the entanglement that 
they can produce. For example, we will determine which 
are the operators Uab that can create maximally entan- 
gled states. While most of our results are concerned with 
two qubits, we will also show that if we allow them to be 
initially (locally) entangled with some ancillas, one can 
obtain more entanglement, at least, for certain measures 
of entanglement. 

In general, an arbitrary unitary operator acting on 
two qubits can be parameterized in terms of 15 coeffi- 
cients (plus a global phase). Thus, to study the maxi- 
mum entanglement which can be produced in terms of 



all these parameters seems a formidable task. However, 
we will show that one can always decompose Uab = 
{Ua <8) UB)Ud{VA <S) Vb), where Ud has a special form 
that only depends on three parameters and the rest are 
local unitary operators. This implies that we can restrict 
ourselves to characterize operators in the form Ud- The 
use of the magic basis introduced in Ref. will also 
considerably simplify our derivations. 

This paper is divided into four sections. In Section 
II we introduce our notation and recall some measures 
of entanglement and their properties. In Section III we 
show that there exists a decomposition of any two-qubit 
gate, which allows us to simplify the problem. In Section 
IV we consider the problem of two qubits. We determine 
how much entanglement can be produced by a general 
two-qubit gate acting on a product state. We also find 
which of such states give rise to that amount of entangle- 
ment. In Section V we discuss the case where we allow 
the qubits to be initially entangled with ancillas. We 
will show that the solution to the problem depends on 
the measure of entanglement we use to quantify it. We 
will also give two examples in which this problem can be 
solved analytically for a particular measure of entangle- 
ment. 



II. DEFINITIONS AND PROPERTIES OF 
ENTANGLEMENT MEASURES 

The purpose of this Section is two fold. On the one 
hand, we give the definitions and notations that will be 
used throughout the whole paper. On the other hand, we 
review some measures of entanglement (for pure states) 
and some of their properties. 



A. Definitions 

We consider two partners, Alice and Bob, who posses 
two quantum systems, A and B, respectively. These sys- 
tems will be composed of one or two qubits each. We will 
express the states of these qubits in terms of the compu- 
tational basis, {| 0), I 1)}. The Hilbert space of system A 
(B) will be denoted by Ha (Hb) respectively. 

Throughout this paper we use capital Greek letters for 
joint states of systems A and B and small letters for 
states describing either system A or system B. We de- 
note by I 5'"'^) a state which is orthogonal to | Vl'), whereas 
I ^f*) denotes the complex conjugate of | ^) in the com- 
putational basis. We will denote the Pauli operators by 
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axjCTy, (Tz and hy a = {gx ,(Ty,az)- If it is not clear which 
system an operator is acting on, we specify it with either 
a sub- or superscript, e.g. a a or a^- 

For two qubits, the Bell basis is defined as follows: 
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= -^(1 00) ± 1 11)), I vl/±) = -^(1 01) ± 1 10)). (1) 

We also make use of the so-called magic basis , which 
is defined in the same way as the Bell-basis except for 
some global phases. We will denote the elements of this 
basis by 



^i) = 1$^ 
$3) = I 



I $2) = -«!*-) 

|$4) =-i|*+) 



(2a) 
(2b) 



The coefficients of a general state in that basis will be typ- 
ically denoted by fik', that is, we write | 5") — J2k ^"=1 'I'a;)- 
In what follows | $) denotes a maximally entangled state. 



B. Measures of entanglement 

We review here some measures of entanglement for 
pure states. In the first part of the paper we are go- 
ing to use the so-called concurrence , C. It is defined 
as 



Ci\^)) = \{^\ay(^ay\^*)\. 
Writing | vj/) in the magic basis we get 



(3) 



(4) 



In the second part of the paper we will use other 
measures of entanglement, which are better expressed 
in terms of the Schmidt coefficients. A pure state | 4*), 
describing the state of two particles, A and B, each of 
dimension m always has a Schmidt decomposition in the 
form: 



VI/) 



^Ck\<l>k)A\tpk)B, 
fe=l 



(5) 



where (0^ |0/) = (-0^ 1-0;) = Ski^kJ = 1, . . . , m. The real 
and positive coefficients Ck , which are the square roots of 
the eigenvalues of the reduced density operator, pA — 
trsd *)(^' I) (or pb = tr^(| *)(*!)) are called Schmidt 
coefficients. We will choose them in increasing order, i.e. 

Cl < C2 . . . < Cm. 

The Entropy of entanglement is defined as follows: 



EE{\^)) = SipA)^-tr[pAl0g,ipA)] 



E 

k=l 



l0g2(Cfc) 



This measure has a well defined meaning: given n copies 
of a state | \1/) then one can produce, using only lo- 
cal operations and classical communication, uEeH ^P)) 



maximally entangled states and vice versa (in the limit 
n — > 00). 

Another useful measure is the Schmidt number [ p^ , 
which we will denote by Es- It is the number of Schmidt 
coefficients which are different than zero minus one. 

There are other set of measures, the so-called Entan- 
glement Monotones, which arise in the context of allowed 
modification of entangled states under local operations 
p3[ . They are defined as 



£;„(!*)) = ^4, 



(7) 



fc=i 



for n = 1, TO — 1. 



We will also use the so-called 2-Entropy (or 2-Renyi 
entropy) of the reduced density operator |14|. It is de- 
fined as 

m 

Sfld *)) ^ Sr{pa) = 1 - tr(p^) = 1 - ^ cl (8) 

fe=i 

where again pA denotes the reduced density operator of 
I *)(^' I and Ck the Schmidt coefficients of | *). In the 
following we will call this measure the Renyi entangle- 
ment. 



C. Properties 

A state describing two qubits contain s two S chmidt co- 
efficients at most, ci,C2, where C2 = y'l — . Thus, its 
entanglement is completely determined by one parame- 
ter, Cl. All the measures of entanglement are monotonic 
functions of each other and, therefore, equivalent. In 
higher dimensions (to-IcvcI systems) though, this is no 
longer true. Let us denote now by | and | v]/') two 
states describing two TO-level systems. Then it might 
happen that for some measure | ^I/) is more entangled 
than I vj/') , whereas for some other measure it is the other 
way around. 

Let us briefly recall some of the properties which have 
to be satisfied by any measure of entanglement, E ||l5|| : 

(a) Monotonicity under local operations: Suppose that 
Alice makes a measurement on her qubit and she 
obtains with probability pk the state ak . Then the 
entanglement cannot increase on average, i.e: 



Eip) > Y,pkE{ak). 



(9) 



(b) Convexity: The entanglement decreases if we discard 
(6) some information, i.e. 



E 



^PkPk 



< 



^PkE{pk 



(10) 
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Now we briefly summarize some useful properties of 
the particular measures of entanglement mentioned in 
the previous subsection. Let us start with the properties 
of the concurrence, C, assuming that we have two qubits: 

(i) Cil'f)) = US fil^ e^^l^fep (fc = 1,...,4). This 

means that a state, written in the magic basis is 
maximally entangled iff its coefficients are real, ex- 
cept for a global phase. 

(ii) C(| *)) = iff I is a product state iff J2k = 0- 

These two properties imply that if | and | ^'^) are 
real in the magic basis (and therefore they are maximally 
entangled) then the state | $) ±i| is a product state. 

Let us also review some properties of the Entropy of 
Entanglement, Ee, the Schmidt number, Es-, the entan- 
glement monotones and the Renyi entanglement, En, 
for arbitrary states of two m-level systems: 

(1) A maximally entangled state, | of two m— level 

systems has m Schmidt coefficients, which are all 
1/^. Thus Ee{\ *)) = log2(m), Es{\ ^)) ^ m - 
1, En{\ 5')) = riEi = n/m, {n — 1, . . . , m — 1) and 
ER{\-^)) = l-l/m. 

(2) A product state can be written as | 5*) = | c^)a| 

thus it has only one Schmidt coefficient, which 
is equal to 1. And so £^£(1*)) = Es{\'^)) = 
= £^i?(|*)) -0 (n= l,...,77i-l). 

III. UNITARY OPERATIONS 

In the next sections we will calculate the maximum 
attainable entanglement produced by two-qubit gates. 
In this Section we consider an arbitrary unitary operator 
Uab acting on two qubits and derive some properties 
which will simplify the problem. 

In Appendix A we show that for any unitary operator 
Uab there exist local unitary operators, Ua, Ub, Va, Vb, 
and a unitary operator Ud such that 

Uab = UA(E>UBUdVA(g>VB, (11) 

where 

Ud = e-''"^^^'''"« (12) 

and d is a diagonal matrix. Here, oa^ denotes the trans- 
pose of OA expressed in the computational basis. We will 
denote the diagonal elements of d by a^, Oy, az- Note 
that any measure of entanglement is not changed by lo- 
cal unitary operators. Thus the entanglement created by 
Uab is the same as the one created by UdVA Vb- And 
so the maximal amount of entanglement which can be 
produced by applying a general unitary, Uab is the same 
as the one created by Ud- This means that we have to 
deal with unitaries which are determined by only 3 pa- 
rameters, (ofa;, ttj,, Qfz), instead of 15 parameters, which 



are used in order to describe a general unitary operator 
acting on two qubits. 

Furthermore, in Appendix B we show that when study- 
ing the maximum entanglement created by a two-qubit 
gate we can restrict ourselves to the case where 

7r/4 >a^>ay>az>0. (13) 

This is due to the fact that the maximal amount of en- 
tanglement created by Ud is symmetric around 7r/4 and 
7r/2-periodic in ax,ay and a^. 

It can be easily shown that the operator Ud is diagonal 
in the magic basis, and therefore we can write 

4 

C/d = ^e-'^'=|*fe)(*fc|- (14) 

The phases are 

Xi = ax ~ ay + az, (15) 
A2 = -ax + ay + az, 
A3 = -ax - ay - az, 
A4 = ax + ay — az- 

IV. TWO QUBITS 

In this section we consider the following scenario. Al- 
ice and Bob have one qubit each. They want to entangle 
them by applying a given unitary operation Uab- Their 
main goal is to find the best separable (pure |l^) in- 
put state that gives as much entanglement as possible. 
According to our previous discussions, we just have to 
find which states \(j))A,\''P}B maximize the concurrence 
of the output state Ud \ 4>) a, \iI')b, where Ud is given in 
(|l^ ) with restrictions (|l3|). We will call these states best 
input states. 

Writing the input and output state in the magic basis 
with the coefficients Wk , respectively, we apply the 
unitary operator Ud and obtain 

J^Mfel ^k) = Ud{\ q^)A\i')B) - J2 ^kc-'^'l^k)- (16) 

k k 

We want to maximize the concurrence of the output 
state, C = IX^fe/^fcl' where we have to make sure that 
the following conditions are satisfied: 

(cl) J2k l/^fcP = 1' which is that the output state is nor- 
malized. Note that, since Ud is unitary this implies 
that the input state is normalized. 

(c2) J2k MfeC^*'^'' = 0. This condition is due to the fact 
that the input state is a product state, which can 
be seen as follows. From Eq.([l^) we see that Wk = 
fikc"^^^, and according to Section IIC (ii) this last 
one is a product state iff the sum of the coefficients 
in the magic basis squared vanishes, which implies 
the above equation. 
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We can determine the maximum of the concurrence of 
the output state under the conditions (cl), (c2) by max- 
imizing and imposing the above conditions in terms 
of Lagrange multiphers, i.e. we maximize 



k,l k 



where 771 is real. We find it convenient to denote 
= Ce'\ 772 = \m\e'\ and = \^ik\e'i'' ■ We 

obtain 



IJ-k^if^'t) 



- my-k + V2f^ke 



Vfc. 



(18) 



Multiplying Eq by /ifc, summing over k and using 
(cl) and (c2) we find that 771 — C^. And so we have, 
assuming that C 7^ 0, 



(19) 



One of the solutions to this equation is fj,k — 0. To find 
the others we write Eq. ([l^) as 



1 _ lMpi(2Afc-7+£) 

c 



c. 



(20) 



Let us distinguish now two cases, namely when 772 is 
zero or not. 

• 772 = 0: From Eq. (|^) it follows that C = 1. Using 
then Eq. it is easy to see that e^''^'^ = e"*''Vfc. 
Thus all the coefficients have the same phase (ex- 
cept for the sign) and therefore the output state is, 
according to the discussions of subsection II C, a 
maximally entangled state. In order to obtain this 
state by applying Ud to a product state the con- 
ditions, (cl) and (c2) still have to be imposed. In 
Appendix C wc show that those conditions can be 
fulfilled iff + ay > 7r/4 and at the same time 
cty + az < 7r/4. There, we also determine the best 
input state. 

• 772 7^ 0: In this case Eq. (|2^) can have at most two 
solutions for a fixed value of I772I/C Thus, in order 
to fulfill ( ^ Vfc at least two of the coefficients have 
to vanish |[l7|]. Let us call the other two fj,k and /i;. 
Then, in order to fulfill conditions (cl) and (c2) 
we have to satisfy l/Xfcp + |^;|2e2''(>'<+«i-(^'= = 
and the normalization condition. Thus |/ife| — 

= I/V2 and the difference between the two 
phases £,k and is A; — — tt/2. With all that it is 
now simple to determine that the largest reachable 
concurrence is 



C = maxfc^; I sin(Afc - A; 



(21) 



Except for global phases the corresponding output 
state is 1/V2(| $fc)+j| $i)e-^'=^^') and the separable 
input state, which leads to this maximum is 



V2 



(22) 



Note that the input state 1/V2(|$fe) - i| $/)) 
(the corresponding output state would then be 
l/v^(|$fc) - i\<^i)e^''-^')) leads to the same 
amount of entanglement. 

Note that in case < 7r/8, we obtain that C = 
sin(a2; + ay), which is directly related to the entangle- 
ment capability of the Hamiltonian of the form adders 
Q . For higher values of ax the result may not be directly 
related to that quantity. 

In summary, in this subsection we have shown that if 
we apply Ud to a separable input state and calculate the 
maximum of the concurrence of the output state, then 
we find: If a^ + ay > 7r/4 and ay + az < 7r/4 then 
this maximum is equal to 1. Otherwise it is given by 
(pl|). In addition we determined the best input state in 
each of those two cases. Note that, since we were dealing 
with two-qubit states we could have taken, according 
to the discussions in section II C, any other measure of 
entanglement to obtain the same result. 



V. USING ANCILLAS 

We analyze now whether and how it would be possi- 
ble to increase the amount of entanglement of the output 
state with the help of auxiliary systems. So, we consider 
the situation where Alice and Bob have two qubits each 
jlSf . Let us denote the auxiliary qubits by A' and B' . 
We allow input states in which Alice and Bob's qubits 
are locally entangled, i.e. of the form | (/>)aa' | V')b-B' ■ 
Then they apply a non-local unitary transformation, 
Uab to the qubits A and B. The question is, then, 
for which \(I))aa' and |V-')-BS' they are able to reach 
max|0)^^,j^)^^,£;(C/^B (g) tA'B'\(t>)AA'A '^)bb'), where 
E denotes some measure of entanglement. In what fol- 
lows we write again simply ( | V') ) instead of | 4>) aa' 
(I 4')bb')- On the other hand, according to Section III, 
we can restrict ourselves to operators Ud of the form (p^). 
For convenience we will call the input state where | 0) and 
I V') are both maximally entangled, local maximally en- 
tangled and the one where both are product states, local 
product state. 

The main difference with the previous section is that 
now the best input states depend on the measure of en- 
tanglement. To illustrate this fact we show in the first 
subsection that for some measures of entanglement the 
best input states are the ones where | (/>) and | ■0) are 
entangled. On the other hand, there are measures of en- 
tanglement for which a local product state is the best 
input state. 
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In the second subsection we show that for some special 
class of Ud (where d has only one non vanishing element) 
the solution to our problem is independent of the mea- 
sure of entanglement. In particular, we show how much 
entanglement can be created in this case and what is the 
best input state. Furthermore, for the class of Uj, where 
all the diagonal elements of d are the same we will de- 
termine the maximum Renyi entanglement as well as the 
best input state according to this measure of entangle- 
ment. 



A. Dependence on the measure of entanglement 

Let us compare the answer to our problem for some of 
the measures of entanglement which we recalled in Sec- 
tion II. According to some numerical examples we have 
that: 

• Schmidt number: Local maximum entangled states 
are always better than local product states. This 
can be easily understood since in the first case the 

maximum value which Es can take is 4. whereas in 
the latter one it can be at most 2. Thus using this 
measure of entanglement the ancillas will in general 
increase the entanglement of the output state. 

• Renyi entanglement: We have checkfxl that for this 
measure the best input states are always either 
local product states or local maximally entangled 
states. In particular, in the next subsection we will 
provide analytical results for some particular cases. 

• Entanglement monotones: We have verified that 
there are unitary operators Ud for which local prod- 
uct states are the best input states, whereas for 
some other values the local maximally entangled 
states lead to the most entangled output state. But 
there also exist some Ud for which neither the prod- 
uct state nor the maximally entangled state are the 
best input states. 

From these examples it becomes clear that it does not 
make much sense to ask for the best input state, if one 
does not specify according to which measure of entangle- 
ment. 



B. Examples 

Before we start with the examples let us make some 

general statement about the input state. It can always 
be written in the Schmidt decomposition as 

I <P)aa' = ca\ 0o)a| 0)a' + sa\ ^o)a\ 1) A' , (23a) 
I 4')bb' = sb\ V'o)i3| 0)b' + Cb\ 4>^)b\ (23b) 

where + = cl + si = 1. This is due to the fact 
that local unitaries applied to A' and B' do not change 



the entanglement (and commute with Ud)- Let us now 
treat two cases in which it is possible to determine the 
best input state for some measures of entanglement. The 
first one should be viewed as a very simple illustration, 
whereas the second one is much more involved. 



1. Example 1 

Let us consider the following simple unitary operator, 

Ud = e'^"^^ = cos(q)1 - i sin(a)ax a-^. (24) 

In this case it is fairly simple to determine the output 
state. It has at most 2 Schmidt coefficients and therefore 
the state can be viewed as a state describing two qubits. 
This implies, as discussed in II C, that all the measures 
of entanglement are equivalent when calculating the op- 
timal states. We take Ee- 

Let us define pi as the density operator whose off- 
diagonal elements are zero, whereas the diagonal ele- 
ments are the same as the one of the state p. Using 
the fact that the von Neumann entropy is convex we have 
that S{pi) > S{p). Apart from that, since the problem is 
symmetric under exchanging system (AA') with [BB'), 
it is easy to verify that the states with ax\(j)) oc !</>"'") 
and CTxIV') oc |f/'"'") lead to the most entangled output 
state. Now, since the states | 0) = 1 1), | -0) = 1 1), as well 
as the states | 0) = | 10) = ! fulfill this condi- 
tion, both, a local product state and a local maximally 
entangled state are the best input states. The maximal 
entropy of entanglement which can be obtained is then 

ma,xEE = — cos(a)^ log2[cos(Q;)^] — sin(a)^ log2[sin(a)^]. 

(25) 



2. Example 2 

Here we determine the best input state, according to 
the Renyi entanglement, corresponding to a unitary of 
the form 

Ud = e-*«(^«+^''+^») = (cos(a)3 -ism{af)l (26) 
-i sin(a) cos(a)e'"(5^ + Sy + S^), 

where Sp = ® {(3 = x,y,z). Here we have used 
[Sx,Sy] = [Sx-Sz] = [Sy,Sz] = 0. lu Appcudlx D we 
show that, according to any measure of entanglement the 
best input state can always be written as 

I <P)aa' = cal 0)^1 0) A' + Sal 1)a| I) A' (27a) 
I V')sB' = Sb\ 0)b| 0)b' + Cb\ 1)b\ (27b) 

wheres^+c^ =s2 + cf = 1. 

Let us start by proving that the best input state is 
either a local product state or a local maximally entan- 
gled state. Calculating the reduced density operator one 
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finds that Pa = Pi ® P2 (tr(pip2) = 0), where pi,p2 
are 2x2- matrices which depend on Sa,Sfc and a. It is 
straightforward to calculate Eii{Ud\ 4') \ tp)) = SbXpa) = 
Snipi) + Sr{p2) and determine its maxima. One finds 
that the either the local product states (s^,Sj = 0,1) 
or the local maximally entangled states = = 1/2 
always lead to a maximum of the Renyi entanglement. 
In the case of a local product state, it is easy to check 
that the best one is | 01) (or equivalently 1 10) @). Let 
us denote by E^'^ {E^^) the Renyi entanglement for 
a local maximally entangled input state (product state 
I 01)) respectively. We obtain 



E'^-^i^a) = ^6^^ - 2cos(4a) - cos(4a)2] 
i?P-^(a) = i[l-cos(4a)2]. 



(28a) 
(28b) 



Comparing those two expressions we find that ao = 
arccos(l/5)/4 w O.IOOtt. So we have that Va < ao the 
local product state is the best input state and otherwise 
the local maximally entangled state leads to the output 
state with the largest Renyi entanglement. In Fig. 1 we 
illustrate this result. 



0.75 




Tt/16 a Ji/8 a 3n/16 nlA 

FIG. 1. Renyi entanglement for the local maximally entan- 
gled input state (1) and for the product state | 01) (2) 

VI. CONCLUSIONS 

We have shown which separable pure two-qubit states 
have to be used in order to create as much entanglement 
as possible by applying a general two-qubit gate. We 
have shown which unitary operators are able to create 
a maximally entangled state. For all the other unitary 
operators we have given the maximal amount of entan- 
glement which can be created by them (Eq.(^l|)). Fur- 
thermore we have shown that by using ancillas one has 
to specify which is the measure of entanglement to be 
maximized. We have given two examples of unitary op- 
erations for which it is possible to determine the maximal 
amount of some particular measure of entanglement. 
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APPENDIX A: DECOMPOSITION OF UNITARY 
OPERATORS 

In this appendix we show that for any unitary op- 
erator Uab acting on two qubits there exist local uni- 
tary operators Ua,Ub,Va,Vb and a non-local unitary 
Ud = e"*'^-^'''^^ {d diagonal) such that 



Uab = Ua® UsUdVA <E) Vb- 



(Al) 



Our prove will be constructive. Let us call a basis consist- 
ing of maximally entangled orthonormal states a maxi- 
mally entangled basis. In what follows the use of the 
subscript k implies that the definition or statement is 
true for A: = 1, . . . , 4, if not stated differently. 

Lemma 1: For any maximally entangled basis {| ^'fc)}, 
there exist phases Cfe and local unitaries Ua, Ub such that 



Ua «) UBe"^" I *A 



(A2) 



Proof: According to the discussion in Subsection IIC 



s'^'H *fc), where I *i 



(i) we can always write 
real in the magic basis. Let us consider two different 
states, I §fe) and | ^i), then 1/V2(| ^k) - i\ ^i)) = I e, /) 
and 1/V2(J §fc) + i\ §;)) - | i, /), where | e), \i) £ Ha 
and I /),!/) G Hb- Note that | e, /) must be orthogo- 
nal to |e, /). This immediately implies that these vec- 
tors must give the Schmidt decomposition of both | ^k,i)- 
Thus we can write 



1 



(A3) 
(A4) 



Using the same arguments for | ^'3,4) 
mine that they can be written as 



it is easy to deter- 



|§3) = ^(e^'|e,/^)-He-V^,/)) 



|*4)=±^(e^1e,n-e-^|e^,/)) 

for some S. In this case, choosing 

UA = \0){e\ + \l){e^\e^\ 



UB = \0){f\ + \l){f' 



(A5a) 
(A5b) 



(A6a) 
(A6b) 
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and the phases Cfc appropriately, we have (|A^). □. 

Note that this first lemma implies that one can go from 
one maximally entangled basis to any other using only 
local unitaries, if one chooses the global phases appropri- 
ately. 

Lemma 2: Given a general unitary operator, then 
there always exist phases ej, and two maximally entangled 
basis {| '^k)} and {| ^'fe)} such that 



(A7) 



Proof: We give a constructive proof. Let us denote 
by {l^'fc)} the eigenstates of U'^U, where denotes 
the transpose of U in the magic basis and e^"*" are the 
corresponding eigenvalues. Note that the eigenvectors of 
the symmetric operator U"^U are orthonormal and real, 
except for global phases. Thus, since we are working in 
the magic basis, they build a maximally entangled basis. 
Now we define I ^'fc) as 



(A8) 



Since the set {| also form an orthonormal basis, 
it remains to prove that its elements are real. In or- 
der to show that let us consider the eigenvalue equation. 



1)1 = 0. Multiplying it by [/*e-"'= we 
0, which is true iff 



get that {e-""'U - e""'U*)\ 
g-iefc[/| ig real □. 

With all that we are now in the position to show that 
any unitary operator can be decomposed into local op- 
erators and Ud as in Eq.(Al). So let us now give the 
procedure to determine the unitary operators that ap- 
pear there. 

1) Calculate the eigensystem of the unitary, symmetric 
operator U^U. Let us denote the eigenvalues by e^*^'' and 
the eigenstates by | '^k)- As proven in Lemma 2 the set 
of those states is a maximally entangled basis. 

2) Choose Va,Vb and the phases ^k, as explained in 
Lemma 1, such that 



3) Calculate 



(A9) 



(AlO) 



Note that according to Lemma 2 the set of those states 
is also a maximally entangled basis. 

4) Choose the eigenvalues of Ud, e''^'' (note that this is 
equivalent to choose the diagonal elements of d) and the 
unitary operators Ua, Ub such that 



C/t (gi uU''-^''^^''+"'M'^k) = I $s 



(All) 



which, according to Lemma 1, is always possible. It is 
simple to check that with these definitions we obtain the 
decomposition (Al). 



APPENDIX B: PERIODICITY AND SYMMETRY 
OF THE MAXIMAL AMOUNT OF 
ENTANGLEMENT 

Let us start out by proving the periodicity of the en- 
tanglement created by Ud- We define d,{d') as a matrix 
whose diagonal elements are ax,ay, ctz (a^ +7r/2, ay,az) 
respectively. It is simple to verify that Ud = —iSxUd'- 
Since Sx is a tensor product of two local unitary opera- 
tors the entanglement created by Ud is the same as the 
one created by Ud'- The same argumentation holds for 
ay and and therefore the amount of entanglement cre- 
ated by Ud is 7r/2 periodic in ax, ay and az- 

To prove the symmetry around tt/A in ax, ay, az of the 
maximal amount of entanglement we use the following 
definition; d, (d') is a matrix whose diagonal elements 
are tt/A + ax,ay,az (7r/4 — ax,ay,az)- It is straight- 
forward to show that Ud = —ia^U^,a^ , where U^, de- 
notes the complex conjugate of Ud' in the standard ba- 
sis. And so we have that E{Ud\'i')) = E{U;i,a^\^)), 
where we used that local unitary operators do not change 
the entanglement. Now, we use that for any measure 
of entanglement, E, i;(| *)) = E{\^*)). This is ob- 
vious, since all the measures are determined by the 
Schmidt coefficients and they are real. Thus, we have 
that E{Ud\ *)) = E{Ud'{(Tx\ *>)*)■ It is clear that the 
maximal amount of entanglement created by Ud is the 
same as the one created by Ud'- Again the same argu- 
mentation holds for the other angles, which proves the 
statement. 



APPENDIX C: TWO-QUBIT GATES WHICH 
CREATE MAXIMALLY ENTANGLED STATES 

We are going to prove here that there exists a nor- 
malized product state | 0)| V') such that | $) — Ud\(f>) \ V") 
is a maximally entangled state iff ax + ay > 7r/4 and 
< 7r/4. 

According to our discussions in Subsection II C and 
Section IV, this is equivalent to fulfilling the conditions 
(cl) and (c2), where = \fik\'^e~'^'^ . Multiplying condi- 
tion (c2) by £-'(^+2^3)^ 

we obtain 



ay + a 



0, 



(Cl) 



where we have defined ai = A{ax ct2 — 'iictx + cez) 

and as — 4:{ay + az)- Note that since tt /A > ax > ay > 
az > 0, we have that 2Tr > ai > a2 > a^ > 0. 
Let us distinguish the following two cases now: 

• ai < TT (or as > n). In this case all the imaginary 
parts appearing in Eq. ( |Cl] ) are positive (negative) 
and therefore the sum can never vanish. 

• ai > TT and a^ < tt. Here the imaginary part of 
|/X4pe'"i is negative, whereas the one of |M2pe'"'' is 
positive and therefore it is always possible to find a 



7 



solution to Eq. (|C1|). In particular we can choose 
/zi = 0. Th en writing the real and imaginary part 
of Eq. (CI) and the normalization condition (cl) 
we simply have to solve: 



sin(a3)|/i2p 
|/i3p + cos(a3)|At2p - 



' sin(Q;i)|/X4p = 
- cos(ai)|/X4p = 



= 1. 



(C2) 
(C3) 
(C4) 



operator u^, 
and cr-'l0n ) 



I CTj^, , where n' is defined as a^, | (/>o) = | 0) 
: — 1 1). Combining all that we have that 



E{UAB\<i>)\^))<E{UAB\{Ca\m 



s,|ll))(sf,|00)+c,|ll))]}, (D£ 



Note that since we have found the solution for the /ifc's, 
it is easy to determine the input state by using the for- 
mula Wk — /ifce*^*-'. 



APPENDIX D: BEST INPUT STATE FOR 
EXAMPLE 2 

Here we prove that the input state which leads to the 
most entangled output state can be written as 



10) =ca|00) + s,|ll) 
|i/^) =Sfc|00)+Cb|ll) 



(Dla) 
(Dlb) 



where + = + = 1. We wiU use that [cr^ ® 
o'ff ,Ud] = 0, where an — a -n. This can be easily verified 
using the commutation relations of the Pauli operators. 

Let us now recall that the input state in system AA' 
can be written as | 0) = Cq| (/)o)a| 0)^' + Sa\ 4>o)a\ 1)a', 
where + s'^ — 1. It is clear that there exists a vector 
n such that cr,-j| (Pq) = | 0o) and (Tfi \ (p^) — — | iPq)- Note 
that I d>) is invariant under cr4 (g) .i.e. 



cr^ (E) I 



(D2) 



Using the fact that Ud commutes with (g) a? and 
with local operators acting on the auxiliary systems to- 
gether with Eq.(p2|) we have that, cr^ (g) af ® ® 
VB'Ud\4>)\'4!) = Ud\ (t))cri ® Vb'I V'), for any unitary op- 
erator Vb'- 

Let us now introduce a new auxiliary system which 
we denote by C. Then, using the convexity (|l^) of any 
measure of entanglement, i5, we have that 

E{Uab\4')aMbb'c)>E{Uab\4>)aMbb'), (D3) 



where 

\i^)BB'C = l/V2(|^)sB' 



0)c + a| 



•Vb' 



\^)bb'\1)c)- 
(D4) 



Now choosing Vb' = <yf and requiring that (g) 
Vb'\iP)bb' — |V')s-B'i which implies that \ip) = 
st\ 0oO) + cb\ we get th&t E{Uab\ (p) aa'\iP) bb') > 
E{UAB\4>)AA'\-tp)BB'), V|0),|'!A), where both, | (^)and 
I ip) are invariant under the operation cr^j (g) tr^. Using the 
same argumentation as before, we can apply the local 
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